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Moments of convex distribution functions 
and completely alternating sequences 

Alexander Gnedin'^ and Jim Pitman*^ 

Mathematisch Instituut and University of California, Berkeley 

Abstract: We solve the moment problem for convex distribution functions on 
[0, 1] in terms of completely alternating sequences. This complements a recent 
solution of this problem by Diaconis and Freedman, and relates this work to 
the Lcvy-Khintchine formula for the Laplace transform of a subordinator, and 
to regenerative composition structures. 

1. Introduction 

It is well known that the distribution function of a probability measure on the 
unit interval [0, 1] is uniquely determined by its sequence of moments 



(1.1) c(n) := / x"-dF{x), n^O,l,... 

J[o,i] 

defined by Lebesgue-Stieltjes integration. A complete characterization of such mo- 
ment sequences was discovered by Hausdorff [15]. To recall his result, for j = 
0, 1, 2, ... , let denote the jth iterate of the difference operator 



Vc(ri) := c(ri) - c(n + 1), n = 0,1,2, 



That is. 



(1.2) V^c(n) =^(-l)M-^)c(n + i), n > 0, j > 0. 



J 
i 

Definition 1.1. A sequence c = (c(0), c(l), . . .) is completely monotone if 
(1.3) VJc(n) > for aU j > 0, n > 0. 

Theorem 1.2 (Hausdorff [I ")]). A sequence c can be represented as the moment 
sequence (1.1) for some probability distribution F on [0, 1] if and only if c{0) = 1 
and c is completely monotone. 
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A sequence c generates a triangular array (c(n, to); < to < n, ti = 0, 1, . . .) by 
the formula 

(1.4) c(n, to) := (^"J V"-"c(to) , < to < n. 

Obviously, c is completely monotone if and only if 

(1.5) c{n,m) > 0, < to < ri. 

For such a sequence c with c(0) = 1, for each n = 0, 1, . . . the numbers (c(n, to), < 
TO < n) form the probability distribution of a random variable Sn with values 
0, 1, . . . ,n: 

(1.6) c(n, to) = P(S'„ = ?7i), < TO < n. 
Hausdorff showed that then there is the convergence in distribution, 

(1.7) lim F{Sn/n < x) ^ F{x), 

n — ^oo 

at all continuity points x of the unique probability distribution function F whose 
Tith moment is c(n). Moreover, 

(1.8) c(n,TO)=f"^/ x"'(l~xV-"'dF(x) , 0<m<n. 

VW J[o,i] 

This probabilistic interpretation of Hausdorff 's theorem can be further developed 
as follows. 

Theorem 1.3 (B. dc Finctti's theorem [!)], §VII.4). Fore completely monotone with 
c(0) = 1, the triangular array (1.4) defines the probability distribution of the random 
number Sn of successes in the first n trials of an infinite exchangeable sequence of 
successes and failures, according to (1.6). For such an exchangeable sequence, Sn/n 
converges almost surely to a random variable X with distribution function F, and 
conditionally given X , the trials are independent with success probability X . 

In a recent study of moment problems [7]. Diaconis and Freedman considered 
the family of probability distributions F on [0, 1] with a nondecreasing density / 
on the semi-open interval [0, 1[ , so that 

F{x) = / f{u)du, < a; < 1. 

J[0,x] 

Note that F{1 — ) < 1 is allowed, in which case the distribution has an atom of 
magnitude 1 — F{1—) at 1. Equivalently, F{x) is a convex function of a; G [0,1], 
with F{0) ~ 0. Diaconis and Freedman showed the following. 

Theorem 1.4 ([7], Theorem 10). A sequence c with c(0) ~ 1 is the moment se- 
quence of a convex probability distribution function F on [0, 1] with F{0) ~ if 
and only if for each fixed n = 1,2,... the sequence c{n,m) is non-negative and 
nondecreasing in to; 

(1.9) c(n, 0) > and c(n, to + 1) — c(n, to) > , < to < n. 
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In view of (1.6) and (1.7), condition (1.9) is a natural discrete analogue of the 
nondccreasing density condition for F. 

In the course of our work on partition structures derived from regenerative ran- 
dom sets [12, 13, 14], we faced a similar problem of characterizing the Laplace 
exponent of a subordinator in terms of its values at positive integers. We noticed 
that this problem was equivalent to the moment problem for convex distributions 
on [0, 1], and that both problems can be reduced to a known integral representation 
of sequences a subject to the following condition, which was studied by Choquet. 
See [1, Definition 6.1, p. 130] for history and terminology. 

Definition 1.5. A sequence a is completely alternating if the sequence — Va is 
completely monotone. That is to say: 

(1.10) V-'a(n)<0, j = l, 2,..., 71 = 0,1,..., 

The fundamental representation of completely alternating sequences is provided 
by the following theorem. 

Theorem 1.6 (Contained in [1], Proposition 6.12). A sequence a is completely 
alternating if and only if there exists a bounded measure v on [0, 1[ and constants 
A >Q and B such that 

f 1 - 

(1.11) a(n) = An + B+ i^(dx) . 

Jlo.i[ 1-2; 

The data (A, B, v) are uniquely determined. 

In Section 2 we show how Theorem 1.6 follows easily from Hausdorff's Theorem 
1.2, and then deduce the following variant of Theorem 1.4. 

Theorem 1.7. A sequence c is the sequence of moments of a convex distribution 
function on [0, 1] with F(0) = if and only if c(0) = 1 and the sequence a defined 
by 

(1.12) a(0) = 0, a(n) = nc(n- 1), n = l,2,... 

is completely alternating. 

Comparing Theorems 1.4 and 1.7, we deduce that condition (1.9) on a sequence 
c with c(0) = 1 is equivalent to the condition that a derived from c via (1.12) 
is completely alternating. We check this directly by algebra in Section 3, thereby 
providing a new proof of Theorem 1.4. We explain in Section 4 how we were first led 
to Theorem 1.7 by consideration of the Levy-Khintchine formula for the Laplace 
transform of a subordinator. Finally, Section 6 relates distributions with higher 
convexity properties to alternating sequences of higher order. 

2. Proofs of Theorems 1.6 and 1.7 

Proof of Theorem 1.6. Let a be a completely alternating sequence with a(0) = 
0. Then — Va is a completely monotone sequence which can be represented by 
Hausdorff's theorem as 



(2.1) 



-Va(n)= / ev{di), n = 0, 1, 
J[oa] 
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for a unique bounded measure v on [0, 1]. Hence, by summing a geometric series 

f 1 - f" 

(2.2) ain)^niy{l}+ ^ v{dO, n = 0, 1,.... 



[oa[ 1 - C 

This is the special case of (1-11) with B — Q. This case, applied after first subtract- 
ing a(0) from aU terms of a completely alternating sequence a, gives the general 
form (1.11) with B = a(0). Conversely, if a is defined by (1.11) or its special case 
(2.2), then (2.1) is obtained by subtraction, hence — Va is the completely monotone 
moment sequence of v. □ 

We note in passing that the linear term nv{l} in (2.5) could be absorbed into the 
integral by extending the integral from [0, 1[ to [0, 1], with evaluation of the inte- 
grand by continuity at ^ = 1. But in this and the similar expressions (1.11) and (2.2) 
we prefer to display this term separately to avoid any possible misunderstanding. 

Our proof of Theorem 1.7 is based on the integral representation of probability 
distributions with convex distribution functions, used also by Diaconis and Freed- 
man. Let 

i^^(x):=^ l(x>0, eG[0,l[, 
which is the distribution function of the uniform distribution on 1] , with density 

Let Fi{x) = l(a: = 1). corresponding to a unit mass at 1. 
Lemma 2.1 ([ ], Lemma 2). The formula 



(2.3) F{x) = / F^{x)v{dO, 

sets up a bijection between convex probability distribution functions F on [0, 1] with 
F{0) = and probability measures v on [0, 1]. This relation between F and v implies 




(2.4) / g{x)dF{x)= / g{x)dF^{x)\v{dO 

"'[0,1] "'[0,1] \ "'[0,1] 

for every non-negative measurable g. 

Proof. A convex distribution function F with F{Q) = has on [0, 1[ a nondccreasing 
density /, a version of which is the right derivative of F. Hence it is clear that a 
unique probability measure v is defined by 

m-l x6[o,i[, 

J[a.x] ^ ^ K 

and v{l} = F{\) — F(l— ), from which (2.4) follows by Fubini's theorem. □ 

Lemma 2.2. Let F be the mixture of F^^ 's with respect to i^id^) for some probability 
measure v on [0, 1], as in Lemma 2.1. Then the moments of F are determined by 

(2.5) n[ x''-^dF{x)=nu{l}+ [ i—^„{d^)^ n=l,2,.... 

J[0,1[ "'[0,l[ 

Let a(0) = and let a{n) be the common value of both sides of (2.5) for n = 1,2,.... 
Then a is completely alternating. Indeed — Va is the completely monotone moment 
sequence ofv, as in (2.1). 
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Proof. By elementary integration, the moments of are given by 

(2.6) n[ x^'-UF^ix) ^ ^-—^ , n=l,2,.... 

J[o,i] 

So (2.5) is the instance of (2.4) with g{x) = nx"~^. The rest is read from the 
previous discussion of (2.2) and (2.1), as formahzed in Theorem 1.6. □ 

Proof of Theorem 1.7. If c is the moment sequence of F^ and a is derived from c via 
(1.12), then a is completely alternating by the previous lemma. Conversely, given 
a sequence c with c(0) = 1, let a be derived from c via (1.12). If a is completely 
alternating, then Theorem 1.6 represents a{n) by the right side of (2.5) for n = 
1,2,.. ., for some unique probability measure v on [0, 1]. Then, by (2.5), the sequence 
of moments of F is the same as if F were the convex distribution function uniquely 
associated with v via (2.4). Finally, F equals this convex distribution function, by 
uniqueness of the solution of the Hausdorff moment problem. □ 

Remark. A subtle point of the above argument is that possibly /jq ^x"^ dF{x) = 
00, in which case the left side of (2.5) has no meaning for n = 0. We insist in any 
case that a(0) = in Theorem 1.7 and Lemma 2.2 by definition. Later discussion in 
Section 4 makes it clear that is the limiting value of A /j^ dF{x) as A | 0. 

But this fact is not relevant to the present argument. 

3. Some algebra 

Here we check directly that condition (1.9) on a sequence c with c(0) = 1 is equiv- 
alent to the condition that a derived from c via (1.12) is completely alternating. 
According to the Leibnitz rule from the calculus of finite differences, the product 
xy of two sequences x and y has successive differences 

(3.1) VJ(x?/)(n) ^ (■^^V^-'x{n + i)V'yin). 

Applied to a{n) := nc{n — 1), with c(— 1) := 0, this gives 

(3.2) V^a(n) = nV^c{n - 1) - jV^'"^c(n). 
A simple computation using (1.4) and (3.2) now shows that 

(3.3) c(n,m+ 1) - c(77.,m) = — I V""™a(m + 1) , < m < n. 

TO + 1 \TO / 

The equivalence of (1.9) and the completely alternating condition on a is now 
evident from (3.2) for n = and (3.3). 

4. Subordinators 

We explain in this section how we were led to formulate Theorem 1.7 by the ap- 
pearance of completely alternating sequences and monotone densities in another 
context: the Levy-Khintchine formula for the Laplace transform of a subordinator. 

Let H denote a probability distribution on [0, oo]. Recall that H is called in- 
finitely divisible if for every n = 2, 3, . . . there exists a sequence of n independent 
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and identically distributed random variables, with values in [0, oo], whose sum is 
distributed according to H. It is well known that H is infinitely divisible with 
-ff[0, cxd[> if and only if H is the distribution of Y{1) for some subordinator 
{Y(t), t > 0), that is a process with stationary independent nonnegative incre- 
ments with values in [0, oo], and Y{0) = 0. Note that this process may jump to the 
terminal value cxo in finite time, but is not identically oo. 

Theorem 4.1 (Levy-Khintchine formula [2]). The Laplace transform of a subor- 
dinator {Y{t), t > 0) is given by the formula 

(4.1) E[e-^^(*)] =e-**(^) 
with Laplace exponent 

(4.2) $(A) = \d+ [ (1 - e-^y)A{dy) 

J]0,oo] 

for some d > and some measure A on ]0, oo], with (d. A) determined uniquely by 
the distribution ofY[t) for any fixed t > 0. 

Remark. The constant d accounts for the continuous drift component of the sub- 
ordinator, while A, called the Levy measure, governs the Poisson rates of jumps of 
various sizes. In particular, the mass A{oo} is the rate at which the process jumps 
to oo. 

Proof. We sketch only the following derivation of (4.1) and (4.2), adapted from [2, 
pp. 5-7], because it highlights the connection with monotone density problems. Let 
(y(t), t > 0) be a subordinator. By independence and stationarity of increments 

E [e-^^(^+*)] ^ E [e-^^(^)] E [e-^^^], 

which implies (4.1) with < ^'(A) < oo by virtue of the well known functional 
equation for the exponential function. Using again the properties of increments, 

$(A) = lim 77z(l - e-*(^)/") = lim mE[l - e^^^fi/")] 

m — 'OO m— i-oo 

= A lim / e^^y mf{Y{l/m)>y)dy , 

where TOP(y(l/m) > y)dy is a measure with nonincrcasing density for each m = 
1,2,.... By continuity properties of the Laplace transform these measures must 
converge, as m oo, to a measure A(y)dy with nonincrcasing density and a possible 
atom d at 0. Thus 

(4.3) $(A) = Ad + A / e-^y'K{y)dy 

J]0,oo[ 

from which (4.2) follows by integration by parts, with measure A defined through 
A]?/, oo] = A(?;) , y > 0, for a right-continuous choice of the monotone density 
K{y). □ 

To relate the Levy-Khintchine formula to the discussion of the previous sections, 
make the change of variables x = e^y in (4.3) to rewrite this equation as 

(4.4) $(A) = Ad + A/ x^^^~K{-\ogx)dx , 

J]o,i[ 
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or again 



(4.5) 



$(A) = A 



10,1] 



where F{x) is the convex distribution function with F{0) = 0, with increasing 
derivative A(— logx) for x € [0, 1[ and with an atom of magnitude d at 1. Formula 
(4.5) shows that the problem of characterizing moment sequences of convex prob- 
ability distribution functions on [0, 1] with F(0) = is identical to the problem of 
characterizing the sequence of evaluations ^'(1), ^'(2), . . . for the Laplace exponent 
$ of a subordinator subject to the normalization condition <&(!) = 1. Since (4.2) 
can be rewritten as 



where i'{dx) is the image of (1 — e ^)A(dz) via x = e ^ , we deduce the following 
corollary of Theorem 1.6. 

Corollary 4.2. There exists a subordinator with Laplace exponent $ subject to a 
prescribed sequence of values 



if and only if this sequence is completely alternating. The value of $(A) is then 
determined for all X > by the Levy-Khintchine formula (4.2), for some uniquely 
determined d and A. 

Theorem 1.7 can be read from Corollary 4.2 and the alternate representation 
(4.5) of the Laplace exponent. This is how we first recognized that the completely 
alternating condition of Theorem 1.7 was equivalent to the condition (1.9) in The- 
orem 1.4. 

Further remarks. According to a classical theorem of Miinz [IS], a sufficient con- 
dition for density in C[0, 1] of the space of functions spanned on a; i— s- x^", n = 
1,2,..., is that A„ — > oo and J2n ^/^n = oo. So the value of $(A) at every real 
A > is determined by the values $(A„) at points of such a sequence (A„). Newton's 
interpolation series [10] allows to explicitly recover 3>(A) from the discrete evalua- 
tions <i>(A„) , n = 1, 2, .... In the same vein as Corollary 4.2, Hausdorff 's theorem 
implies that a sequence Vl'(O) = 1, ^'(1), ^'(2), ... is the sequence of evaluations of 



for some distribution of X on [0, oo] if and only if this sequence is completely 
monotone. In that case the distribution of X, and hence ^'(A) for all A > 0, 
is uniquely determined. As observed by Feller [s], consideration of the sequences 
"^{0) = 1, *(e), *(2e), ... for e = 2"'' and fc = 0, 1, 2, . . . , leads to Bernstein's theo- 
rem that a function ^' is the Laplace transform of some probability distribution on 
[0, oo[ if and only if * is infinitely differentiable with ^{0) = 1 and (-l)"v]/(") > Q 
for every n, where ^E*'") is the nth derivative of 

Similarly, a function $ is the Laplace exponent of some subordinator if and only 
if <t> is infinitely differentiable with <i>(0) and (-1)"(()(") < for every n. See [9] 
for further discussion, and [1] for the theory of completely monotone and completely 
alternating functions defined on a semigroup instead of the positive integers or the 
positive halfline. 



(4.6) 




<I>(0) = 0, $(1), $(2),... 



*(A) =E(e--^^) 
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5. Regenerative composition structures 

Wc sketch in this section a probabilistic interpretation of completely alternating 
sequences, based on our recent work [12, 13, 14] on regenerative compositions and 
their associated partition structures. We note in passing that developments and 
applications of de Finetti's theorem have played an important role in a number of 
studies of measure- valued processes and associated particle systems [4, 5, O]. The 
general idea is that some random allocation or splitting of a mass continuum can 
be described in terms of a simpler combinatorial model obtained by independent 
random sampling of n points in the continuum. Probabilities in the combinatorial 
model are typically represented as moments of random variables of interest in the 
continuum model. The combinatorial models are consistent in a natural sense as 
the sample size n varies, and the continuum model is recovered as a law of large 
numbers limit. 

The simplest illustration of this idea is provided by de Finetti's theorem for 
sequences of zeros and ones (Theorem 1.3). Suppose the unit interval is split by a 
point X into two interval components [0, X] and ]X, 1]. Let Ui, U2, ... be a sequence 
of independent uniform [0, 1] variables independent of X, and let 

n 
i=l 

Associating each Ui with a 'ball', and [0, X] and \X, 1] with two 'boxes', Sn describes 
an allocation of n balls in two boxes, with Sn balls in the left-hand box, and n — Sn 
balls in the right-hand box. As n varies, these allocations are sampling consistent, 
meaning that if a ball is picked uniformly at random and deleted from the nth 
random allocation, the result is distributed like the (n — l)th random allocation. To 
paraphrase de Finetti's theorem: every sampling consistent sequence of distributions 
for the allocation of balls Sn can be realized via this scheme directed by some 
random variable X S [0, 1], and X is recovered as the limit of Sn/n as n ^ 00. 

A straightforward generalization of this model provides an interpretation of de 
Finetti's theorem for exchangeable sequences with a finite or countably infinite 
number of possible values. Let 

= Xo < Xi < X2 < • • • < 1 

with lim„ Xn = 1 be the cumulative random sums 

k 

i=i 

associated with some random discrete distribution (Pi, P2, ■ • •)■ Let !7i, [/2, • ■ • be a 
sequence of independent uniform [0, 1] variables independent of the X^s, and set 

n 

Nn^j := ^ l(X,_i <U,< X,), J = 1, 2, ... . 

i=l 

Then the random sequence of counts 



{Nn,l,Nn,2, ■ ■ .) with J2 ^".J = " 
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can be interpreted as a random allocation of n balls into a sequence of boxes labeled 
by J = 1, 2, . . .. These allocations are sampling consistent in an obvious sense, and 
according to de Finetti's theorem every sampling consistent sequence of allocations 
can be represented this way, with Pj the limiting frequency of balls in box j. 

In some applications of random discrete distributions, the labeling of atoms of 
the distribution is of little or no importance. If all that matters is the relative sizes 
of these atoms, and perhaps some ordering of these atoms (not necessarily a simple 
indexing by positive integers), the natural combinatorial object is an allocation 
of n balls into some finite number of non-empty boxes, which might be cither 
ordered or unordered. For n a positive integer, a composition of n is a sequence of 
positive integers with sum n. A partition of n is a non-increasing composition of n. 
The terms of a composition may be called its parts. A composition structure is a 
sequence (C„, n = 1,2,...) of random compositions of integers which is sampling 
consistent, that is: 

If n identical balls are distributed into an ordered series of boxes according to 
(C„), then a distributional copy of C„_i is obtained by discarding one of the 
balls picked uniformly at random, independently of {€„), and then deleting 
an empty box in case one is created. 

This is a variation of Kingman's notion of a partition structure [17], which is a 
sequence of random partitions of integers {Pn) subject to the same consistency 
condition, except that after discarding a ball and deleting an empty box if necessary, 
the boxes are permuted to obtain a partition of n — 1. Kingman [17] established 
a one-to-one correspondence between partition structures (Pn) and distributions 
for a sequence of nonnegative random variables Vi , V2 , . . . with Vi > V2 > . . . and 
X^i Vi < 1. In Kingman's paintbox representation, the random partition P„ of n is 
constructed as follows from (14) and a sequence of independent random variables 
Ui with uniform distribution on [0, 1], where (Ui) and (Vfc) are independent: Pn is 
defined to be the sequence of ranked sizes of blocks of the partition of [n] generated 
by a random equivalence relation ^ on positive integers, with i ~ j if and only if 
either i = j or both Ui and Uj fall in Ik for some k, where the Ik are some disjoint 
random sub-intervals of [0, 1] of lengths Vk- See also [-(}] and papers cited there for 
further background. 

Gnedin [11] gave a similar representation of composition structures, using a ran- 
dom closed TZ C [0, 1] to separate points of a uniform sample into clusters. Given 
TZ, define an interval partition of [0, 1] comprised of gaps, that is open interval com- 
ponents of [0, 1] \ TZ, and of individual points of TZ. A random ordered partition of 
the set [n] := {1, . . . ,n} is constructed from TZ and independent uniform sample 
points Ui, . . . ,Un by grouping the indices of sample points which fall in the same 
gap, and letting the points which hit TZ to be singletons. A random composition 
Cn of n is then constructed as the sequence of block sizes in this partition of [n] , 
ordering the blocks from left to right, according to the location of the corresponding 
sample points in [0,1]. Gnedin showed that every composition structure (C„) can 
be so represented. As in Kingman's representation of partition structures, TZ can 
be interpreted as an asymptotic shape of C„, provided C„ is properly encoded as an 
element of the metric space of closed subsets of [0, 1] with the Hausdorff distance 
function. 

An interesting class of composition structures (C„) is obtained by supposing that 
TZ is the closure of {1 — exp(— F(t)),< > 0} for some subordinator {Y{t),t > 0). 
As shown in [12], these composition structures are characterized by the following 

regenerative property: 



Moments of convex distribution functions 



39 



For all n > m > 1, given that the first part of C„ is m, the remaining 
composition of n — m is distributed like Cn-m- 

Let q(n, m) denote the probability that the first part of (C„) is of size m. Then for 
each composition A ~ (Ai, . . . , Ac) of ri 

£ 

(5.1) F{Cn = X) = l[q{A,,X,), 

i=i 

where Aj = Aj + . . . + Xg. We showed in [12] that this formula together with a 
recursion implied by sampling consistency forces 

(5.2) ,(n,™) = ^HL_^ ) 

for some unique sequence $(n), n = 0, 1, 2, . . . with $(0) = and $(1) = 1. Non- 
negativity of the matrix q shows that $ is completely alternating. So Corollary 4.2 
now provides the integral representation 



(5.3) <^{n) = nd+ (1 - (1 - a;)")£>(dx) . 

J]o,i] 

for i>{dx) the image of the Levy measure A{dy) of the subordinator via the trans- 
formation from y to 1 — exp(— y), and hence 

(5.4) g(7i,m) = ndl(m= 1)+ f") / a;'"(l - a;)"-™z>(dx). 



m 



0,1] 



6. Higher convexity 



A distribution function F on [0, 1] is said to be (fc + \)-convex if the derivative 
p(k-i) exists and is a convex function on [0,1[ (see [3, 16, 19] and references 
therein for this and other generalized concepts of convexity and their applications 
in statistics). The moment problem for this class of distributions can be analyzed 
in the same way as for the convex distributions we considered above (which were 
2-convex, with k = 1). For a (fc + l)-convex F, Taylor's formula with remainder 
becomes 

(6.1) F{x)=P{x)+[ i^^^(dO, xe[0,l[, 

i[o,x] (1-4) 

where v is a, bounded measure on [0, 1] and P is a polynomial of degree at most 
k — \. We suppose further that P = 0. In this case the analogue of (2.5) becomes 

{n + lpc{n) ^A{n + if^ + fc! / 

J[o.i[ 

where nJ'^ ~ n{n + 1) • • • (n + fc — 1) and A > 0. A principal role is delegated to the 
k-associated sequence 

a(0) = . . . = flfc-i = 0, a/j(n) = n*^^ c(n — fc) for n > k, 

which is k-alternating sequence meaning that (— V)'^afc is completely monotone. A 
generalization of Proposition 1.7 emerges. 



fc-1 

E 

i=0 



{x - ly 



v{dx) 

[x-iY 
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Proposition 6.1. For each A; = 1, 2, ... , a distribution function on [0, 1] is (fc + 1)- 

convex with F^^\Q+) = for j = 0, . . . ,k — 1 if and only if c(0) = 1 and the 
k-associated sequence Ok is k- alternating. 

The side condition in tlic proposition kills the polynomial part P in (6.1), a 
condition which was guaranteed by F{0) = in the case A; = 1. If the coefhcients 
of P are nonnegative, it is still possible to give a characterization in terms of a 
fc-associated sequence alone, whose initial k terms are no longer zeroes, rather need 
to be defined in a nontrivial way by Newton's interpolation. 
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